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The scattering of ultrasound from bubbles of ~ l/xm radius, such 
as used in contrast enhancers for ultrasound diagnostics, is studied. We 
show that sound scattering and "active" emission of sound from oscil- 
lating bubbles are not contradictory, but are just two different aspects 
derived from the same physics. Treating the bubble as a nonlinear oscil- 
lator, we arrive at general formulas for scattering and absorption cross 
sections. We show that several well-known formulas are recovered in 
the linear limit of this ansatz. In the case of strongly nonlinear oscil- 
lations, however, the cross sections can be larger than those for linear 
response by several orders of magnitude. The major part of the inci- 
dent sound energy is then converted into emitted sound, unlike what 
happens in the linear case, where the absorption cross sections exceed 
the scattering cross sections. 

PACS numbers: 87.59.Mt, 43.80.+p, 43.25.+y, 43.35.+d 

I. INTRODUCTION 

In recent years, a new kind of contrast agent for use in ultrasound diagnostics has 
been developed: suspensions of gas filled microbubbles, i.e., bubbles of at most a few 
micrometers in diameter The enhancement of image brightness and contrast is 
mostly due to the well-known fact that microbubbles of this size can have tremen- 
dous scattering cross sections for the incident diagnostic ultrasound (with frequencies 
around 1 — lOMHz) A number of different theoretical approaches are present in 
literature, resulting in formulas for the scattering and absorption cross sections that 
do not always seem compatible (see e.g. |jl]J^-§). It is one of the main goals of this 
work to show that all of these results can be understood as special cases of a unified 
approach which treats the bubble as a (generally nonlinear) volume oscillator. The 
different appearance in literature is only due to the use of different formalisms (e.g. 
full-fledged scattering theory with partial wave decomposition |^ or linear oscillator 
theory |Q) and the validity of different limiting cases (e.g. neglecting surface tension 
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The other main goal of our present study is to go beyond previous work [0, 
where monochromatic driving was treated, and employ pulsed driving within our 
formalism, in order to achieve a more realistic modeling of the situation in diagnostic 
applications. The analysis of the emitted sound from the bubbles in terms of intensity 
and spectral distribution is of obvious importance to assess the signal quality of 
ultrasonography with bubble contrast agents. While a large number of publications 
have dealt with the sound emission of linearly oscillating bubbles (see e.g. and 
references therein), and the weakly nonlinear case was treated by Prosperetti [^^, a 
systematic study of cross sections for strongly nonlinear response to pulsed driving 
has, to our knowledge, not yet been performed. It is this fully nonlinear case of 
bubble oscillations which is encountered in ultrasound diagnostics as driving pressure 
amplitudes of up to ~ 10 atm are common, whereas strong nonlinearities occur already 
for amplitudes > 1 atm. 

We present our general, unifying approach in Section ^ and show in Section |T| 
how it translates to a multitude of linear limiting cases, where analytical results can 
be obtained if monochromatic driving is assumed. Still in the linear limit, we then 
introduce pulsed driving in |^ and finally present results for the full nonlinear case in 
Section 0. A more detailed discussion, especially of spectral properties of the emitted 
sound for nonlinear response, will be postponed to another publication [^]. Section [VI| 
presents conclusions. 



II. SOUND EMISSION AND ABSORPTION 

A. Emitted sound pressure 

We want to evaluate the pressure of emitted sound Ps{r,t) from a body capable 
of volume oscillations (i.e., a bubble) driven by an incident pressure wave P{t). First, 
we notice that there are two contributions to Ps{r, t): (i) the active emission of sound, 
caused by the change of volume of the body, and (ii) the passive contribution due to 
the mere presence of the (maybe non-oscillating) body in the incident field. Thus, we 
can write 

Ps{r,t) = P:{r,t) + P!{r,t). (1) 

For a body of volume Vi,{t) which is much smaller than the wavelength of the 
incident sound, the pressure Pf(r, t) of actively emitted sound is given to leading 
order by P,p!0| 



'^^ Attt dt^ ' ^ ^ 

Here, pi is the liquid density and we have adopted a spherical coordinate system 
with radius r. There is no angular dependence because we have assumed the long 
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wavelength limit (S wave scattering). For a spherical bubble with radius dynamics 
R(t), (H) easily translates to 

P;(r, t) = i2R{tf + R{t)m) = -q:it) , (3) 

where we have introduced = rPg{r,t), thus separating the trivial (geometrical) 
1/r spatial dependence of the radiation from the time-dependent part. 

The second term Pf (r, t) in (|1]) is the pressure of passively emitted sound due to 
the mere presence of the body. It arises from the perturbation of the density field 
in water: if the body was not present, the incident sound wave could induce density 
changes in the volume of liquid occupied by the body. As this cannot happen, the 
situation is equivalent to having an oscillator that compensates the effects of the 
sound wave; this oscillator has a (virtual) volume change 

dV dpi/dt 
at pi 

with the (now constant) volume Vq of the body, see 0. For a sound wave P{t), we 
have dP/dpi = cf, where the derivative is taken at constant entropy. Then, in analogy 
to (g), replacing VJ, by V and inserting Vq = inR^/S, we get the leading order pressure 
change due to the passive reaction to the sound wave 

P!{r,t) = ^,R'{t)P{t) = lqm, (5) 

where q is the speed of sound in the liquid. Note again that we are in the limit 
R{t) <^ A, where A is the wavelength of incident sound. Therefore, the driving 
pressure P(t) can be treated as spatially uniform, i.e., the pressure experienced by 
the bubble does not vary over its size. In analogy to (^, we have defined the quantity 



B. Bubble oscillation 

To evaluate (^) or (|^), we need a formula for R{t). The oscillatory behavior of the 
radius P of a gas bubble in a liquid is well described by the Rayleigh-Plesset equation 
TT| , p!0[| . Many variants of this equation have been presented (see e.g. |]12|-|1^); the 



following form has proved to be robust and close to experiment even in situations of 
massively nonlinear bubble behavior, as e.g. in sonoluminescence experiments [0: 



RR +Ir' = - {p{R, t) - Pit) - Po) + —^AR. t) - - ^ . (6) 
2 pi PiCidt R piR 

We have introduced the liquid viscosity ui and surface tension C, here, as well as the 
ambient pressure Pq taken to be 1 atm in this work. All liquid parameters assume 



3 



the values of water, except for z//, which is multiphed by three to mimic the viscosity 
of blood [^]. p{R,t) stands for the gas pressure inside the bubble and is modeled 
by a polytropic process equation with van der Waals hard core and parameters for 
air. The polytropic exponent is taken to be one, because bubbles of the sizes we treat 
here (~ l/zm) are smaller than the thermal diffusion length on the time scales of the 



oscillation, and can therefore be regarded as isothermal [21,22 



C. Scattering cross sections 

In general, the scattering cross section age is related to the incident intensity 
(energy/area/time) Imc and the scattered power Wsc (energy/time) via 

Wsc = (Jsc line , (7) 

and has thus the dimensions of determined by the incident pressure 

wave -P(t), namely (see, e.g., [^, §65) 

Iinc=—{P\t))t (8) 
PlCl 

and Wsc follows from Ps{r,t) through 

Att Att 
Wsc = —{r'P^{r,t)), = —{qmt. (9) 
PlCl piCi 

Here, {■)t denotes a time mean; we have exploited the virial theorem to express 
intensities and powers solely in terms of pressures (the "potential" part of the energy 
of the wave). The integral 

Att /"'" 

Esc = — ql{t)dt (10) 



PlCl JO 

gives the total scattered energy over a time span r (e.g. the duration of an incident 
sound pulse). 

From these definitions, it is clear that the general formula for the scattering cross 
section is 

Att Att 

^ wwy-^^^^' ^ ww}^'^'-^^^ ^ ■ ^^^^ 

We notice that the active and passive parts interfere. It is therefore, strictly speaking, 
wrong to divide cTsc into an active and a passive part, as it is sometimes done in 
literature. Nevertheless, it will be shown below that in the case of diagnostic bubbles 
the total scattering cross section is almost exclusively due to the active part of 
(0), while the passive contribution can safely be neglected. 
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There are other possible contributions to Ps{r,t), e.g. direct or indirect results of 
bubble shape oscillations. These are not analyzed here because perfect sphericity is 
assumed. Also, there is a contribution if the bubble (or the body in general) can be 
translated as a whole by the incident sound. This term is orthogonal to those treated 
above (i.e., the interference terms vanish) and results in a well-defined additional 
scattering cross section (cf. P] §78 or ||l|) with a characteristic angular dependence. 
As it is (for the case of a gas bubble in a liquid) always of the same order of magnitude 
as the passive contribution, and is therefore equally negligible, we do not treat it here 
in detail. Finally, if two or more bubbles come close to each other, they will change the 
emitted sound field either by direct secondary scattering or indirectly by modifying 
their modes of oscillation (e.g. via secondary Bjerknes forces p3|p^ ). We will not try 
to incorporate these effects, but restrict ourselves to the case of a single bubble. 



D. Absorption cross section 

We have so far regarded the scattering cross section as an indicator of the acoustic 
energy deflected by the scatterer. Likewise, the absorption cross section a^^^ stands 
for the energy loss induced by the viscous term in the RP equation (|^). In this 
case the energy is directly converted into heat and will, in general, be of little use 
to the experimenter, whereas the scattered sound can be detected much more easily. 
Nevertheless, a^^^ is important in order to assess the energy balance of the scattering 
process. The liquid viscosity exerts a stress 

4:iyipiR , . 

Pvis = ^ (12) 

over the bubble surface of size A = AttR^ (this can be read off directly from the RP 
equation). As the bubble wall moves with velocity the dissipated power is 

ly;,, = p,i,AR = IGnuipiRR'' . (13) 
The absorption cross section is determined in analogy to o"sc via 

W^dl = <.s/mc. (14) 



III. SMALL DRIVING: THE LINEAR CASE 

The formulas given in the previous sections can be applied to all cases of bubble 
motion and subsequent sound emission. For small driving, they reduce to the linear 
case and almost all quantities can be calculated analytically. First of all, the RP 
equation can be linearized: we set R(t) = Ro{l + x(t)) {Rq is the radius of the 
undriven bubble under ambient conditions) and get 
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2 P{t) 
X + 270; + UJqX = "2 



(15) 



with the viscous damping constant 7 



^ and the bubble eigenfrequency 

-ft,. 



+ 



4C 



(16) 



This latter quantity consists of two distinct terms, the first due to the gas pressure 
p{R,t), and the second governed by surface tension (. For p{R,t), a polytropic ideal 
gas formula was chosen. While it is advisable to employ a more elaborate formula 
(e.g. a van der Waals gas with hard core) if the bubble oscillation is violent, the ideal 
gas is a very good approximation in the linear limit. The polytropic exponent Hg 
measures the gas compressibility {p{R,t) oc Pg^)', it is 1 for isothermal behavior and 
equal to the adiabatic exponent for adiabatic behavior of an ideal gas (e.g. 7/5 for 
air). In deriving Eq. (|15]), we note that effects of thermal damping are neglected (as 
they have not been present in (^, either), and also radiation damping (caused in 
(^) by the term involving q) is not present. Both of these damping effects are very 
small in the parameter range of our interest; thus, 7 only contains viscous damping 
contributions. 

To get explicit solutions, let us assume monochromatic driving P{t) = Pa cos uJdt 
in (plSf). Pa is the driving pressure amplitude, = 271 the angular driving frequency. 
This yields 



x{t) = e cos{ujdt + 6) 



(17) 



with the amplitude 



and the phase shift 



6 = arctan 



27t^d 



UJj - UJn 



(19) 



Obviously, in the monochromatic case /«, 



11) we get 



^. With R{t) from x(t), (D, (D, and 



V2ujj cos(a;rft + S) V2 2^ 2 



(20) 



Note that a^^ does not depend on the absolute size of P^, as long as the oscillation 
stays linear. Fortunately, in all cases of interest to bubble ultrasound diagnostics, i.e.. 
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FIG. 1. Numerical computation (solid line) and scaling behaviors of the normalized scat- 
tering cross section for bubbles driven with a monochromatic 50kHz sound wave at small 
drive amplitudes (linear case) . We chose a representation similar to that of Nishi |5[ , dis- 
playing asc divided by the geometrical cross section ttRq on a double logarithmic plot. To 
the right of the resonance peak (for large Rq), asc/iT^Ro) relaxes to 4 (dot-dashed line). 
In the opposite limit of small Rq, it scales oc Rq (dashed) for an interval, but becomes 
steeper again for very small Rq, reaching an asymptote oc -Rq (long dashed). These scal- 
ing laws correspond to the Rq and Rq behaviors of the unnormalized asc, respectively (cf. 
Eqs. (|1),(|1)). 



for frequencies in the MHz range and Rq between a few tenths of a micron and a few 
microns, the second (passive) term of this formula is negligible compared to the first 
(active) part. We will show in detail below that the dominance of the active emission 
is so overwhelming that the relative errors in neglecting the passive contribution never 
exceed 0.25% in the parameter range of this study. Dropping the passive term, we get 
the considerably simplified expression (reproducing results presented e.g. in Leighton 
§ in Section 4.1) 

This formula can also be translated into the S wave scattering limit of the complete 
scattering analysis by Nishi 0. In equation (47) of that work, there are additional 
contributions to 7 due to thermal and radiation damping, as well as modifications 
involving the stiffness of the bubble. In certain limiting cases, (|2l| ) gets further 
simplified. Moreover, a rich variety of scaling behaviors, especially in Rq, can be 
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found. The richness (as compared to "ordinary" hnear oscillators) is due to the fact 
that in (|l^), the damping, driving frequency and dimensionless driving all depend on 
Rq. In Fig. |I| cr^*" is presented for a 50kHz-driven bubble, a frequency chosen such 
that most of the limiting cases treated in the following subsections can be illustrated 
by the graph. Note that all cross section graphs in this work were computed with the 
full numerical formulas presented in Section 0. 



A. Small bubbles: dominance of gas pressure or surface tension 

For small bubbles (such that ujq 3> ujd) the denominator of the active contribution 
in (|20|) simplifies to Uq. It is crucial here to verify which of the two contributions to 
ujq in (|TB[) is dominant. If surface tension is very small (and i?o is not extremely tiny), 
uq will be governed by the first term, i.e., the one due to the internal gas pressure. 
With the polytropic ansatz, we can use = HgPo/pg for the speed of sound in the 
gas Cg and from ( pOD the limiting behavior for small Rq is then 

with the compressibilities Ki = l/{picf) for the bubble interior {i = g) and exterior 
{i = I). This reproduces the equations presented e.g. by Landau |^ §78, deJong |^ or 
Meerbaum in Nanda/Schlief |l| (apart from the translational contribution mentioned 
in [11 (J| above). Note that in the case of gas bubbles in a liquid. Kg will be much 
greater than Ki, so that to very good approximation 



and the dominance of the active scattering is again confirmed. In fact, as for typical 
materials pi ~ lOOOpg and q ~ 5cg, the active emission can surmount what would 
be expected from a purely passive scatterer by more than eight orders of magnitude! 
This is why oscillating bubbles are so much superior to (completely passive) hard 
spheres in terms of ultrasound scattering capability. In Fig. |I| there is indeed a region 
where the oc Rq behavior of cr^*" is observed. Note, however, that the validity of ( p2D 
and ( P3| ) for bubbles is limited, especially as surface tension is explicitly neglected. 

If surface tension is taken into account and the parameters for water are inserted 
into (P^), we see that for Rq < 0.96pm the surface tension term dominates ujq, so 



that — > AC,/ piRq. As a consequence, the active scattering cross section will acquire 
a different limit, namely 

^ 4^2 • (24) 
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The steeper oc Rq behavior takes over in Fig. |^ for very small bubbles, as expected. 
As the passive contributions scale like Rq, there has to be a (small) radius where the 
latter become dominant. This radius is easily calculated to be 



= TT^ , (25) 



which, for the material constants of the water/air (blood/air) system, is ~ 4.4-10 ^^m, 
i.e., on a subatomic scale which is not described by the physics discussed here. 



B. Linear resonance 

When Rq is such that uq ~ Ud, the bubble is near resonance, i.e., its oscillation 
amplitude e and thus a^*" become much larger than for neighboring Rq. E.g., for the 
case of air bubbles in blood driven with a;^ = 27r ■ 3 MHz, it reaches values of about 15 
times the geometrical cross section ttRq (cf. Fig. ^ below), i.e., the normalized cross 
section cr'*J'/(7ri?Q) is larger than for any other radius. From (pID, it is easy to get 



_/ jDres\6, ,2 



where i?Q*^* is determined by the condition ujq{Rq^'^) = uj^- The height of the resonance 
peak is thus mainly determined by the strength of damping, i.e., the viscosity of the 
liquid. In water with its three times lower viscosity, the maximum cross section is 
therefore 9 times higher than in blood. The (Lorentz shaped) resonance peak is the 
hallmark of the classical pictures of hnear bubble scattering cross sections (see [|l],H,||). 
We will see that this shape changes considerably in the nonlinear case. 



C. Large bubbles: the "soft sphere" limit 

For large bubbles (but still obeying R <^ X), ujd^ ojq and 

a'::^ At^RI (27) 

i.e., the scattering cross section becomes four times the geometrical cross section, 
which is also reproduced in Fig. |I|. Equation ( P?]) describes the behavior for all Rq if 
the bubble interior is arbitrarily compressible while surface tension is absent (C = 0). 
This can be seen if we imagine — > (implying Cg ^ and Kg —>■ oo) and thus 
Uo 0: the resonance peak is shifted to very small Rq and every larger ambient 
radius results in the "soft sphere limit" cross section (|27|). 

In reality, for large bubbles there are of course still the interference contributions 
with the passive term oc Rq. But the influence of Pf can only outweigh the active 



9 



emission for Rq > -^A (which is again easy to prove); this is aheady in a region 
where the assumption Rq <ti X is dubious. In fact, for diagnostic ultrasound with 
Ud = 271 ■ 3MHz, A ~ 500/im, resulting in a critical radius for passive influence 
of Rq > 140/im. Thus, there is at most a very small region of transition until the 
scattering cross section has to be described by completely different formulas (Mie 
scattering theory), and the bubbles in this region are of no interest to ultrasound 
diagnostics because of their size. 



D. Other limits 

There are a number of limiting cases of (|2^) and (^) that go beyond the parameter 
space of this study; they are given here for completeness: 

• In the limit of Rq ^ A, the cross section relaxes to the limit age = 27ri?Q [^. 
This is only relevant for mm-sized or even larger bubbles. 

• If we apply our formulas to solid spheres in a liquid rather than bubbles, i.e., 
we make the interior of the "bubble" much more stiff than the liquid, the main 
contribution in (^) is the passive scattering, as Kg/Ki <^ 1. The resulting 
cross section is the "hard sphere" Rayleigh limit 

In analogy to the case of the "soft sphere" above, the resonance peak is shifted 
to arbitrarily large radii for an arbitrarily stiff bubble, so that here ( 128|) is valid 
for every finite Rq. 



E. The regime of ultrasound diagnostics 

It is important to note that for the typical /im size bubbles of ultrasound diag- 
nostics applications, the relevant limiting cases of the linear formulas pO|), (pl|) are 
(|^),(|^), and (1^). The case of cr^*" oc Rq given by (^2]) does not occur, because 
the resonance radii Rq^'^ are - for MHz driven bubbles - in the range of the crossover 
ambient radius Rq ~ 0.96/im between gas pressure and surface tension dominated 
resonances. Thus, the oscillations of bubbles with Rq < i?™* are all surface tension 
dominated, and never show the Rq behavior. This is obvious in Fig. |^, which shows 
the linear cross sections for 3MHz driven bubbles in blood. In contrast to Fig. |l|, 
there is a direct transition from the resonance peak to a curve a^*" oc Rq. 
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FIG. 2. Numerical computation (solid lines) and scaling behaviors of the normalized 
linear scattering (thick lines) and absorption (thin lines) cross sections for bubbles driven 
with a monochromatic 3MHz sound wave (diagnostic frequency). The resonance features 
as well as the large Rq limit of 4 (dot-dashed) are again present in asc/iT^Ro)- Of the 
small Rq limits, only the surface tension governed cr^*" oc Rq remains (long dashed). The 
corresponding scaling laws for (y^abT I ^'^^o) (dashed) for large and oc Rq (long 

dashed) for small Rq. 



The linear analysis is naturally extended to absorption cross sections; the analyt- 



Here, the limiting cases for small bubbles show (JabT °^ -^o small Rq for gas 
pressure dominated ujq and cr^hT ^ -^o surface tension dominated case (see 

Fig. ||. Thus, 2> erf,*" for small bubbles, i.e., most of the energy incident on the 

bubble is converted into heat through viscous damping, whereas only a small part is 
available for sound emission. The already surprisingly high scattering cross section 
for microbubbles is outnumbered considerably by the absorption cross section. The 
same is true in resonance, where 



F. Linear absorption cross sections 
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can be calculated for = 3MHz and the material constants of blood, in very good 
agreement with numerical computation (cf. Fig. y). 

For large Ro above the resonance radius we have (JatT -^o^S a'*" grows like 
i?Q, we have in this range cr^^" <^ cr^^c^ and the major part of the energy goes into 
sound emission. 



IV. PULSED DRIVING 

A. The incident pulse 

One assumption that allowed us to give simple analytical formulas was the 
monochromaticity of the driving. In reality, bubbles in ultrasound diagnostics are 
driven by the signal of a diagnostic transducer, which is almost always pulsed, and 
quite often (e.g. in pulse wave doppler mode) the pulses are only a few wavelengths 
long. 

In closer agreement with experimental reality, we therefore model the pressure 
pulse P(t) as 

P(t) = PaCosMt-t,))expi -A^t-Q^\ , (31) 

centered around tc with relative width h. We choose h = 1/3 here. After Fourier 
transform, the frequency space representation is 



{iU-iUd) \ , / {uJ + LUdY 



(32) 



This spectrum is (almost) Gaussian in shape; the corresponding power oc |P(c<j)p 
decays to 1/e^ of its maximum value within a distance of i/icurf around Ud (cf. Fig. |^). 
Thus, the reciprocal 1/h oi the relative width is an approximate measure for the 
spatial extension of the pulse (in wavelengths). In our example, the pulse is about 3 
wavelengths long and corresponds to the shortest pulses routinely available in medical 
applications of diagnostic ultrasound. 



B. Effects on the cross sections 



As we are still in the linear limit, the pulsed driving can be understood as a 
superposition of monochromatic waves of frequency uj with amplitudes proportional 
to the Fourier components P{oj) in (^2]). The response will be therefore obtained 
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FIG. 3. Modeled ultrasound diagnostics driving pulse. The upper part of the figure shows 
the time series P{t) {t normalized to T = 1/fd), the lower part the normalized Fourier 
spectrum of the acoustic power given by P{t). It is centered around the mean frequency 
fd = 'i MHz and decays to 1/e^ of its maximum value within hfd = fd/3 = 1 MHz above and 
below fd- This corresponds to the effective spatial extent of the pulse of 3 wavelengths. 

from a convolution of the monochromatic response discussed above and the spectrum 
P(^). 

Figure | shows the cross sections (Jsc(-Ro) and cr^ft^(i?o) for monochromatic and 
pulsed driving in comparison. The figure is a blow-up of the region near -Rq*^*, which 
is the only area where there are marked differences between the responses to the 
different drivings. As it may have been expected, the response curves for polychro- 
matic driving are broader, because unlike the monochromatic case, where the single 
frequency Ud corresponds to a single, well-defined Rq'^^, there are a number of Rq for 
which the bubbles react strongly to the most intense frequency components in P{uj)- 
Correspondingly, a bubble driven at resonance by a monochromatic wave is a more 
effective scatterer than for pulsed driving, so that the maximum height of the cross 
section curves is smaller in the latter case. 

Moreover, with pulsed driving there is a slight, but significant shift in the maxi- 
mum of both cross sections towards larger Rq. This is due to the asymmetry in the 
response curve around the resonance radius that is apparent already in the monochro- 
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matic case (cf. Fig. on increasing Rq from -Rq*^**, the cross sections do not drop 
as rapidly as when Rq is decreased from the maximum of the curves. Therefore, 
the contributions in P{uj) with u shghtly smaller than Ud (which correspond to 
Ro{u!) > RQ^^ioJd)) will have a larger effect on asc and cr^^^ than the contributions with 
uj slightly smaller than uod-, even if they are represented with the same weight in P{oo). 
Consequently, upon convolution of the spectrum, the cross sections are increased for 
the larger Rq and the maximum of the curves is shifted towards Rq > R^'^'^lujd). 
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FIG. 4. Scattering (thick lines) and absorption (thin lines) cross sections for bubbles 
driven by a monochromatic 3 MHz sound wave (dashed) and a diagnostic pulse according 



to (31) with center frequency 3 MHz (solid). The response to the polychromatic driving 
displays a broader, lower resonance peak, which is shifted towards larger Rq. Further away 
from resonance, the differences to the monochromatic case soon become negligible. 

If the width of the pulse is varied, the peak shift and peak broadening vary cor- 
respondingly: thus, if the pulse becomes very short (shorter than one wavelength), 
the resonance structure is almost completely lost. Using longer and longer wave- 
trains (smaller h), on the other hand, leads to successively sharper resonance peaks, 
while the repetitive oscillations of the bubble allow for the occurrence of subharmonic 
resonances, provided that the bubbles remain shape stable (cf. 0). 



V. RESULTS OF FULL NONLINEAR COMPUTATIONS 

Let us now come back to the general formulas presented in Section |T| and drop the 
assumption of linear response. We will find that the features of the nonlinear case can 
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be quite different from the linear results presented above. Not surprisingly, the cross 
sections now depend on the driving pressure amphtude Pa, which is not the case for 
linear driving. It should be noted that, in order to compare these results with exper- 
imental measurements, they have to be convoluted with the bubble size distribution, 
as in virtually all experiments the bubble suspensions are not monodisperse. 



A. Scattering cross sections 



With P{t) given by pT|), line can be computed, so that with (|T|), (H) and 
results. Let us first convince ourselves that the passive part of (P is tiny in the 
nonlinear case, too. In Fig. |^ we present the relative errors made in determining 
scattering cross sections from P^{r, t) (or qg{t)) alone rather than from the full inter- 
ference formula (|TT]). It is easy to see that these errors are nowhere greater than about 
0.25%, and in most cases much smaller. The huge advantage of (actively emitting) 
bubbles as scatterers compared to stiff solid bodies of similar size is here once more 
strikingly demonstrated. For computational simplicity, the cross sections presented 
in the following figures were computed using only the active emission part. 
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FIG. 5. Relative error of scattering cross section computations with and without passive 
emission taken into account. The numerical calculations were performed using the com- 
plete formula ( pT] ) with or without the passive contribution and driving with an ultrasound 
diagnostics pulse according to (|3T[). The curves correspond to driving pressure amplitudes 
Pa of 1 — 10 atm, in steps of 1 atm. The error is nowhere greater than about 0.25%. 



Fig. ^ shows the normalized cross sections for various Rq G {O.l/im, 5.0/im} and 
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Pa = 1 — lOatm (in 1 atm intervals). These values represent typical bubble sizes 
found in diagnostic bubble suspensions (often the size distribution shows a peak 
around i?o ~ Ifim and typical sound pressure amplitudes in the focus of clinical 
ultrasound devices. From linear analysis, we know what the cross sections for Pq — >■ 
look like. This curve is readily reproduced for small Pa- While the curve for latm is 
still essentially unaltered compared to the linear case, there are dramatic changes in 
(^sc{Rq) with higher Pa- E.g., the resonance radius i?™* undergoes a shift away from 
its linear value to smaller Rq, i.e., in the opposite direction of the shift induced by 
the polychromatic pulse. This shift is explained by the nonlinearity of the oscillator 
(cf. I^Bj): when expanding the RP equation up to nonlinearities of third order, one 
obtains for the nonlinear eigenfrequency 



It is easy to see that with these values we always have ui^ < ujq, i.e., the bubble is 
an oscillator with a soft potential. With stronger nonlinearity, becomes smaller 
for a given Ro; in order to be in resonance, we require Ud = uj^^ with constant Ud- 
Therefore, Rq must be decreased in order to increase c^o beyond Ud to ensure that the 
resonance condition Ud = uJq^ is again fulfilled. Thus, the nonlinear resonance radius 
j^res,NL gj^aller than its linear counterpart. Note also that the resonance structure 
is blurred especially for high driving; this is in contrast to the case of monochromatic 
driving [|], where well-defined resonance radii can be recognized up to the highest Pa- 

A most striking feature of Fig. ^ is the tremendously enlarged scattering cross 
section, especially at radii smaller than the linear resonance radius (i.e., in the region 
of nonlinear resonance). Sometimes, age is greater by several orders of magnitude 
compared to the linear case (note the logarithmic scale). Only for Rq ^ j^res,im 
the shape of the curve virtually unaltered. Thus, the scattering cross sections for 
sma// bubbles (below ^ Rq = Ifim) are severely underestimated hj the linear theory. 

The reason for this effect is found in the typical R{t) radius dynamics of nonlinearly 
driven bubbles: when Pa is sufficiently high, they undergo violent collapsing instead 
of just a smooth oscillation. This leads to high velocities and extraordinarily high 
accelerations of the bubble wall, exceeding lO^g in the solutions of @ for the highest 
Pa of our study (even larger accelerations are known in other contexts of bubble 
dynamics, see [|l9l). These values of R and especially R, which far exceed those 
expected for linear response, lead to very high active emission pressures (cf. Eq. (||)) 
and thus to a huge age- 




(33) 



with the amplitude a of the oscillation and 




(34) 
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FIG. 6. Total scattering cross sections of bubbles driven by ultrasound pulses. Prom 
bottom to top: Pq = 1 — 10 atm in steps of 1 atm. The linear profile of age is drastically 
changed for high Pa, the resonance radius shifts to smaller Rq due to nonlinearities and the 
cross sections are much larger around the resonance radius than in the linear case (up to 
3-4 orders of magnitude). 



B. Absorption cross sections 



In Fig. the viscous absorption cross section cr^^^ computed from ([13|) and (|14D 
is shown for the same parameter combinations as for asc in Fig. ^ As expected, 
the resonance peak undergoes similar nonlinear shifting. However, although a^bg also 
shows a tendency to grow for higher Pa, it is obvious that its dependence on driving 
pressure amplitude is nowhere near as dramatic as in the case of the scattering cross 
section. This is because in the computation of W^^^ via ([T3|) , only the velocity R 
is present, but not the acceleration, which is responsible for the tremendous sound 
emission, as argued in the previous section. As usual, the radiative processes are 
connected with acceleration, while the dissipative processes are governed by velocity. 
For high Pa we can therefore conclude that - regardless of the size of the bubble - 
most of the incident sound energy is scattered again into sound, and not dissipated 
via viscous forces. 
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FIG. 7. Total viscous absorption cross sections of bubbles driven by ultrasound pulses. 
From bottom to top: Pa = 1 — W atm in steps of 1 atm. The change in shape with growing 
Pa is not as drastic as for the scattering cross sections, cr^^^ is considerably smaller than 
(Tsc for the highest Pa, while it is much larger in the linear case of small driving. 

C. Secondary absorption 



It is tempting to conclude that for the highest Pa we have the optimal situation 
for detection of scattered sound. This is not necessarily so, however, because of the 
spectrum of the emitted acoustic radiation. We show in a separate paper that the 
most strongly driven bubbles emit sound in a spectrum of immense band width, with a 
large portion of the energy in the ultra high frequency part (GHz). These frequencies 
are readily absorbed (on a length scale of cm or less) by water or other media (blood, 
tissue) frequently encountered in diagnostic applications (cf. e.g. p7| , p8| ). With this 
absorption, most of the sound energy is converted into heat. Fig. ^ shows the example 
of the effectively detected scattered sound after the pulses have traveled through a 
5 cm layer of water. The absorption properties of water are well known: the energy 
contained in every Fourier component (frequency / = uj/2tt) of the emitted sound 
signal decays over a distance r like 



^,(r,/) = E,(0,/)exp(-«^/V) 



(35) 



with the known characteristic absorption coefficient of water ~ 1.5 ■ 10 ^^s^ jm 



29| . Thus, the highest frequency components experience the strongest damping, and 



it is especially the high-power, high-frequency emission of scatterers at small Rq which 
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is radically diminished by this process, although the cross sections do not drop to or 
near their linear values. The implications for diagnostics, such as possible risks of the 
heat deposition connected with the absorption, are dealt with in separate work 0. 




FIG. 8. Effective scattering cross sections of bubbles driven by ultrasound pulses and 
separated from the detector by a water layer of 5cm width, which acts to damp out high 
frequency sound according to (|35|). Prom bottom to top: Pa = 1 — lOatm in steps of 1 atm. 
Tlie growth for large Pa is much diminished here as compared to Fig. ^. The high frequency 
components of the sound emitted by the small bubbles at high Pa are absorbed in water 
and thus the major part of the emitted sound energy is lost to the detector. 



VI. SUMMARY AND CONCLUSIONS 

We have presented a unified view of gas bubble sound scattering and absorption. 
Starting from the description of the bubble as a nonlinear oscillator, general formulas 
have been developed that can be shown to transform, in certain limiting cases, into 
a number of apparently discrepant formulas found in literature, whether they be 
derived from oscillator theory 0], scattering theory P] or the theory of solid body 
sound scattering It is important to state that there is no antagonism between 
"scattering" and "active sound emission" . Both terms refer to the sound that can be 
detected after incident sound has hit a gas bubble. The total scattering cross section 
does contain a contribution of passive sound emission, but this contribution is tiny and 
can always be neglected for situations in ultrasound diagnostics. Bubbles driven in the 
nonlinear Tange show a much larger scattering cross section than in the linear case and 
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are more "actively" emitting sound in the sense that the incident energy is primarily 
converted into sound and not into viscous heating, as is the case for linear driving. 
The most prodigiously radiating bubbles, at small radii and large driving amplitude, 
emit such high frequency sound that it is again absorbed by the surrounding medium 
and leads to secondary heating. The effective yield of scattered sound energy is 
therefore much smaller than what would be expected without taking absorption into 
account. Still, the ultrasound scattering capabilities of bubbles and therefore their 
effectiveness as a contrast agent in diagnostics are considerably enhanced due to the 
nonlinearity of the oscillations. 
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